Spin and Majorana polarization in topological superconducting wires 
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We study a one-dimensional wire with strong Rashba and Dresselhaus spin-orbit coupling (SOC), 
which supports Majorana fermions when subject to a Zeeman magnetic field and in proximity of 
a superconductor. Using both analytical and numerical techniques we calculate the electronic spin 
texture of the Majorana end states. We find that the spin polarization of these states depends 
on the relative magnitude of the Rashba and Dresselhaus SOC components. Moreover, we define 
and calculate a local "Majorana polarization" and "Majorana density" and argue that they can 
be used as order parameters to characterize the topological transition between the trivial system 
and the system exhibiting Majorana bound modes. We find that the local "Majorana polarization" 
is correlated to the transverse spin polarization, and we propose to test the presence of Majorana 
fermions in a ID system by a spin-polarized density of states measurement. 
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Introduction-MdiioT diiidi fermions have been attracting 
a lot of interest recently in the light of the discovery 
of new materials with topological-insulator properties^. 
These atypical fermionic particles have been predicted 
long time ago by E. Majorana as real solutions of the 
Dirac equation^. Many condensed-matter systems such 
as PfafRan states in fractional quantum Hall (FQH) 
systems^, chiral p-wave superconductor^ (like stron- 
tium ruthenato^), nodal superconductors under cer- 
tain conditions^, ultracold fermionic atoms with laser- 
field-generated spin-orbit intearctions^, the surface of 
3D strong topological insulator^, as well as semi- 
conductor/superconductor heterostructures^ have been 
proposed as platforms supporting Majorana fermions. 
Among the possible heterostructures, one-dimensional 
(ID) systems with a strong spin-orbit coupling such as 
InAs and InSb wires^^, subject to Zeeman magnetic 
field and in proximity of a superconductor (SC), can ex- 
hibit Majorana fermions at their extremitie^i^. Var- 
ious proposals have been made to detect the Majo- 
rana states including interferometry^^, noisei^ and spec- 
troscopy measurements^^, etc. However, while a direct 
confirmation of the existence of the Majorana states 
would constitute an important step for fields such as 
quantum computation^^, they have not been detected ex- 
perimentally so far. 

Majorana modes for two dimensional spin-triplet topo- 
logical superconductors have been shown to exhibit an 
Ising-like spin density that may allow their detection via 
coupling to a magnetic impurity^^. Along similar lines, 
we propose a method to detect the Majorana states in 
ID topological semiconducting wire spectroscopically, us- 
ing spin-polarized scanning tunneling microscopy (STM). 
We generalize the model in Refs. Ellill to include both 
Rashba and Dresselhaus spin-orbit interactions, and we 
show that the resulting Majorana bound states exhibit a 
characteristic spin texture. In particular we find that the 
component of the spin polarization in the transverse spin- 
plane (orthogonal to the direction of the magnetic field) 



is non-zero solely in the topological phase, and its orien- 
tation is determined by the relative weight of the SOC 
components. Moreover, we introduce a "Majorana pseu- 
dospin" local order parameter and define two new local 
quantities denoted Majorana polarization and Majorana 
density which quantify locally the Majorana character of 
a state. We show that the transverse spin polarization is 
related to the Majorana polarization and we propose that 
its measurement via spin-polarized STM will allow one to 
directly visualize the Majorana fermionic states, and thus 
to test the topological character of a one-dimensional sys- 
tem. 

Model- We consider a semiconducting wire oriented 
along the x-direction, and in proximity to a s-wave su- 
perconductor. Due to bulk inversion asymmetry, semi- 
conducting wires can exhibit along with the Rashba SO 
interaction analyzed in Refs. [l2lll8[ a Dresselhaus SO in- 
teraction of the same order of mas; nitudei^ (- 0.1 eV 
A). The Bogoliubov-de Gennes (BdG) Hamiltonian for 
the infinite wire with both types of SOC can be written 
as 



H 



H = y— - /i + apay + (Spa^jTz + VzCTz - Ar^. (1) 

cr's and r's are the usual Pauli matrices acting respec- 
tively in the spin and particle-hole spaces. The chemical 
potential is denoted by /i, Vz is the Zeeman field, A is 
the induced superconducting pairing and a {j3) charac- 
terize the strengths of the Rashba (Dresselhaus) SOC 
components. The presence of the Dresselhaus term only 
trivially modifies the spectrum for the translationally in- 
variant system^ ^ 



(2) 



with ^ = p'^ /2m — ji. A careful analysis of this model 
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shows that the conditions for the existence of the topolog- 
ical phase supporting Majorana fermions are unaffected 
by the Dresselhaus SO interactions, V'^ > + /i^. It is 
interesting to note however that Majorana bound states 
can exist even in the absence of the Rashba term, when 
Dresselhaus SO interactions are present. Most impor- 
tantly, the spin texture of the Majorana states is in- 
fluenced by the presence of the Dresselhaus term. To 
support this claim we present an analytical study of 
the wavefunctions corresponding to the Majorana bound 
states, and we complement it by a numerical study of the 
corresponding lattice model. 

Analytical solution-It has been shown that Majorana 
bound states can arise at the interface between trivial 
and topological regions of a one dimensional wire, for 
example by considering a position dependent chemical 
potent iaii^i^. Thus, by choosing a chemical potential 
/ii for X e [0,L], such that /i^ < — A^, and a 
A^o > ~ outside this interval, one obtains a finite- 
size topological region inside a topologically trivial phase. 
The chemical potentials are chosen such that the p = 
gap, A — Y^K^^^/?, is much smaller than the super- 
conducting gap A, which allows one to obtain analytical 
solutions to the problem by linearizing the Hamiltonian 
in p. We assume that L 1, such that the problem can 
be solved independently at the two ends. Thus, ignoring 
the finite-size effects, the condition to have zero-energy 
solutions bound at the two interfaces yields the allowed 

values for the momenta, ^/oP^^^W^f = A ± ^J'^z ~ j ' 
j G {0, 1}. Consequently, the Majorana solution at the 
left boundary iIj{x ~ 0) can be written as 



, tan^i, , , .-^ 
tan 00 



(1 



tan 01. , . ™ 
— ^ U2 /io e^°" 
tan 00 



X > 0, 



X < 0. 
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Similarly, the Majorana solution at the right boundary 
'0(x ~ L) is 



/^U3(/ii)e" 



-k^ (x — L) 



(l + ^)u3(Mo)e-'^o-(-^) 
tan 00 

(1 - ^)u4Mo)e-^(^-)^ 
tan 00 



X < L 



X > L 



(4) 



with e^^^- = 1/a/2(v'1 + j^j/V^ + i^/l - Iij/V^), and the 
Majorana eigenvectors are given by 

ui (/ij ) = (cos 0j e*^ , — sin 0j , sin 0j e*^ , cos 0j ) , 

U2 (/ij ) = (cos e*^ , sin 0j , — sin 0j e*^ , cos 0j ) , 

U3 (/ij ) = — (cos e^'^ , sin 0j , sin 0^ e*'^ , — cos 0j ) , (5) 

U4^{fij)^ = ( — cos 0j e*'^ , sin 0j , sin 0j e*'^ , cos 0j ) . 



Here e*^ = {a -\- i^^j \J o? -h allows one to define a 
two-dimensional spin-orbit vector in the transverse plane, 
= (cos "d^ sin "S). Note that the obtained wavefunctions 
are indeed Majorana fermions respecting the reality con- 
dition through the phase choice (i9+7r)/2 for the complex 
coefficient k. The magnitude of k is determined from the 
normalization c ondition s of the wave fu nctions and is of 
the order of ( ^A^^T^ - A/ VaM^)^/^ 

The electronic local spin polarization s(x) of a given 
four-component (two spinx electron/hole) Nambu state 
can be calculated by evaluating the expectation val- 
ues Si{x) = {i)\(Ji ^^^^ \^) (the To + Tz insures that we 
only take into account the electronic, and not the hole 
degrees of freedom). This prescription yields for the Ma- 
jorana fermionic states at x = and x = L\ 



s(0) = (— sin(20i) cosi^, sin(20i) sini9, cos(20i)) 

s(L) = ^(sin(20i)cosi^, -sin(20i)sini9,cos(20i)). 

^ (6) 

The above results show that the spin z-components are 
equal at the ends of the wire, while the transverse spin 
polarization is equal in magnitude and opposite. Its di- 
rection is fixed solely by the relative weight of the Rashba 
and Dresselhaus SOC components, Sy/sx = — It 
would be interesting to see whether this results also holds 
in presence of interaction n^^i^^ . 

Definitions- A general wave function written in Nambu 
basis described in Eq. ([1]) as {u^^ u^^v^^ v^) can be recast 
in a Majorana basis (7it^ 72t^ 72;), where 

(715,725) = -^(^^e-^^/2+V^,e^^/^^(V^^e-^^/2-V^,e^^^^ 

(7) 

and is a phase characterizing the particular choice of 
Majorana basis. We can focus simply on the spin up 
part of the wave function and define the Majorana po- 
larization P^f the difference between the probability 
to have a 71^ Majorana and the probability to have a 



72^ Majorana, P^f = 



7it 



P^2t- Note that Pxi/2 but 



not Pm, have been also introduced in Ref. Il8|. For 
the generalized {(^ not fixed) Majorana basis, it reads 
= -2Re[u^^v*e-''^]. The phase Lp defines a "Ma- 



jorana polarization axis" such that P^f can be inter- 
preted as a vector (we can denote this vector as "Ma- 
jorana pseudopsin"). As such it is decomposed on the 
axis defined by ^ into Pm^{^ = 0) = —2Re[u^v^] and 
Pl^iip = Tr/2) = -2Im[u^v;]. 

The same procedure can be repeated for the spin down 
component. This leads to the definition of the full Majo- 
rana polarization, Pm^ = P]^. +^m. with i G {x^y}. We 
should note that the absolute value of the polarization 

vector Vm = \JPm^ + -> which we denote Majorana 

density, is not dependent on the choice of {x^y) axes in 
the Majorana space. For further reference the Majorana 
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polarization and density are 

Pm, = 2Re[uivl - u^v^], = 2Im[^i;v| - u^v^], 

Vm = 2\u^vl-u^v^\. (8) 

To emphasize the utihty of such definitions the follow- 
ing properties should be noted. For a wave function con- 
taining only electronic (or hole) degrees of freedom, the 
Major ana polarization is always zero. It equally vanishes 
for a conventional s-wave superconductor where the wave 
function has the symmetry u^^^^^^ = A non-zero 

value for the Majorana density is a necessary, but not 
sufficient condition to have Majorana fermions, and iden- 
tifies the low energy regime in which the model yields an 
effective p-wave type superconductivity. At zero energy 
this quantity is maximal and predicts Majorana bound 
states at the edges. Inspecting the Majorana polariza- 
tion of these states allows one to unambiguously identify 
that the two Majorana fermions are of different type. 
The most important property is that such definitions al- 
low one to explore locally (on-site in the discretized sys- 
tem) the structure of the wave function, and its Majorana 
character. 

For the above analytical solutions, the Majorana po- 
larization vectors Pm = (^m^ , ) yield 

Pm(0) = -Pm{L) = -|A^p(cos^,sini9cos(2(^i)) (9) 

The Majorana polarization vectors are opposite for the 
two Majorana wave functions at the two ends. Besides, 
one can identify a relation between the spin polarization 
vector and the Majorana polarization. When /i. A, Vz 
are fixed, Pm^ is proportional to Sx, while Pjviy is pro- 
portional to Sy. Thus, when only Rashba/Dresselhaus 
SOC is present, the total transverse spin polarization is 
proportional to the Majorana polarization, with a pro- 
portionality constant which depends on the chemical po- 
tential potential and the applied Zeeman field. When 
both components of the SOC are present, the Majorana 
polarization and the transverse spin polarization vectors 
are no longer collinear but the correlations between these 
two quantities are still qualitatively preserved. 

Numerical model-Foi the numerical study, we consider 
a tight-binding formulation of BdG Hamiltonian in Eq. 

o 



^|(t + iaay + il3(Jx)Tz^j+i + h.( 



(10) 



t _ 



(c]^, c]^, Cj^, -c^-^). The sum 



with the Nambu basis 

is performed over TV = 100 sites in the system. We work 
in units of t = 1 and we consider the lattice constant / = 1 
and h = 1. Numerical simulations are done for typical 
values of the parameters Vz = 0.4, A = 0.3, /i = 0. 

By exact diagonalization we have access to the local 
density of states, and the local spin-polarized density of 



states along the and z directions. For example the 
local (site n) electronic i-spin polarization density is given 

by 



AN 



To + Tz 



l*^f^> (11) 



where Ej is the j*^ eigenvalue of H and \^n ) is 
the site n component of the j*^ eigenvector, |^n^) = 
(ul^^^ul^l^^nl^vl^^). We can thus use the definitions in 
Eq. (jHj) to calculate the local Majorana polarization and 
density. For example, the local (n site) Majorana x- 
polarization can be written as 

PM.,n = E^(^ - ^.)2Re[^.2<f - ulilv^f]. (12) 

Numerically, we implement the S functions as very thin 
Gaussians of width ^ O.OOOIHvf/I' 

Numerical results-As expected, the exact diagonaliza- 
tion of H recovers the Majorana bound states at the ends 
of the wire. In Fig. ^ we present the x and z compo- 
nents of the spin polarization, as well as the Majorana 
polarization when the only considered SOC is Rashba. 
The symmetrical situation, with only the Dresselhaus 
component of the SOC present, is analyzed in the Supple- 
mentary Material (SM). We can see that the zero-energy 
Majorana wavefunctions are extended over a small num- 
ber of edge sites, while exhibiting strongly damped spa- 
tial oscillations. The transverse electronic spin polariza- 
tion is opposite at the two ends of the wire; when only 
the Rashba term is present only the x-component of the 
transverse spin is non-zero. The z-spin polarization is the 
same at the two ends of the wire. While Eq. (|6|) predicts 
that the z-polarization vanishes for /ii = 0, this is not 
the case in the numerical calculation. This can be un- 
derstood by the presence of an effective fi due to the ne- 
glected kinetic term, which to leading order, contributes 
to (p^) ^ 0((A — 14)^), which creates a negative effec- 
tive potential as in the numerical results. Moreover, this 
effective chemical potential is responsible for the spatial 
(quickly damped) oscillations of the spin polarization ob- 
served numerically. Allthough these oscillations are not 
captured by the continuum limit calculations, one can 
however check that the ratio Sy^i/sx^i depends only on 
the spin-orbit couplings in agreement with Eq. (j6j). 

The numerical results for the Majorana polarization 
presented in Fig. ([T]) also follow closely Eq. (|9]). Thus, the 
values of the Majorana polarization are always opposite 
at the two ends of the wire, and only Pm^ is non-zero 
when only the Rashba component is present. Also, as 
expected, Pm^ is proportional in this case to the x-spin 
polarization. Moreover, while we do not present this here 
in detail, when both the Rashba and Dresselhaus SOC 
components are present, the peculiar dependence of the 
Majorana polarization in Eq. ([9j) on the cos(20i) term 
indicates that the Majorana polarization vector should 
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FIG. 1: The spin polarization along the z and x directions, 
and the Majorana polarization Pm^ 5 cis a function of energy 
and position, A = 0.3, Vz = 0.4, a = 0.2, /3 = and /i = 0. 



exhibit a spatial precession in Majorana space; we have 
verified numerically that this is indeed the case. 

We examine if Majorana polarization is a good or- 
der parameter to characterize the topological transition. 
This is done by varying one of the parameters (A, Vz, fi) 
to drive the system in a trivial phase. In Fig. (j2j), we vary 
Vz, and indeed we see that the system becomes trivially 
gapped (no Majorana bound states) for Vz < A. The 
inset describes the dependence of the half-wire integral 
of the Majorana polarization for one of the lowest-energy 
states as a function of Vz (an integral of 0.5 is equivalent 
to a "full" Majorana state). The Majorana polarization 
decreases smoothly to zero below the critical value of Vz . 
The same phenomenon can be observed in the second 
panel, where we plot the spatial distribution of the Ma- 
jorana polarization as a function of Vz. We have noted 
that the transition becomes sharper when increasing the 
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FIG. 2: First panel: lowest-energy eigenvalues and the 
half-wire integral of the Majorana polarization (inset) as a 
function of Vz. Second panel: Majorana polarization of the 
lowest-energy state as a function of position and Vz. 
Parameters: A — 0.3, = 0, /3 = 0, and a — 0.2 



size of the system. The same qualitative features are 
obtained when A and ji are varied across the topological 
transition (the dependence on /i is presented in SM). This 
shows that the Majorana polarization (and density) is a 
good local order parameter for the topological transition 
at V^ = A^ + /i^. This suggests that the Majorana po- 
larization can be used to investigate disordered wires^, 
and indeed, as shown in SM, in the presence of disorder it 
exhibits interesting features such as a weak polarization 
of the low energy bulk states^. Moreover, the spin and 
Majorana polarization exhibit similar spatial structure, 
even in the presence of disorder. 

Conclusion-To summarize, we have found that the Ma- 
jorana end states are oppositely spin-polarized in the 
transverse spin-plane, and the direction of polarization 
depends on the relative weight of the Rashba and Dressel- 
haus SOC in the wire. Moreover, we have proposed a new 
wave- function-based measure of the Majorana character 
of a system, which we denote Majorana polarization. We 
have seen that this quantity is related to the electronic 
spin polarization and we have proposed to test the Majo- 
rana character of a ID system using spin-polarized STM 
measurements. While the density of states measurements 
can only give information about the existence of a lo- 
calized state at a given energy, without telling anything 
about its Majorana character, such a spin-polarized mea- 
surement can make the difference between a Majorana 
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excitation or a non-topological localized state. 
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FIG. 3: In the first panel, the lowest-energy eigenvalues and 

the half-wire Majorana polarization integral (inset) are 
plotted as a function of /i. In the second panel the Majorana 
polarization of the lowest-energy state is plotted as a 
function of position and /x. The parameters considered are 
A = 0.3, = 0.4, a = 0.2, and /3 = 0. 
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Supplementary Material 

The Majorana polarization is a good order parameter to 
characterize the topological transition. This can be seen 
as well when varying the chemical potential fi in Fig. (j3j). 
We analyze also the effects of disorder on the Majorana 
polarization. For example, we consider an on-site po- 
tential Wj which takes random values with a uniform 
distribution in the interval [— ^, A thorough analysis 
of various types of disorder (phase slips, correlated disor- 
der, disorder in the order parameter A and the hopping 
parameter etc.) will be presented elsewhere^^. For a 
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disorder strength of w = 0.5, when the Dresselhaus SOC 
P is absent, the Major ana polarization Pm^ stih behaves 
quahtatively similar to the x-axis spin polarization, as 
described in Fig. (j4]) (^/-polarizations remain zero). Note 
that, for this disorder strength, Majorana states can still 
be observed on the edges. Moreover the bulk states near 
the gap become slightly Majorana- and spin-polarized 
(however they are not full Majorana states). 



Finally, we consider the presence of a Dresselhaus 
SOC, when the Rashba SOC is absent. We find that 
Majorana fermions are present, as illustrated in Fig. (|5]). 
The local density of states reveals zero energy edge 
states, and a plot of the Pm^ identifies them as Majo- 
rana states. Here the transverse spin polarization has 
only a non-zero component along the ?/-axis, that, as 
expected, is proportional to Pm^ • 
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FIG. 4: The dependence of the Majorana polarization and 
x-spin-density on energy and position for a weak-disorder 
w = 0.5. A = 0.3, Vz = 0.4, = 0, a = 0.2, and /3 = 




FIG. 5: Density of states and Majorana polarization Pm^ at 
A 0.3, Vz 0.4, = 0, a = 0, and /3 = 0.2. 



